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SPI MEASUREMENTS OF THE DIFFUSE GALACTIC HARD X-RAY CONTINUUM 
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ABSTRACT 



INTEGRAL Spectrometer SPI data from the first 
year of the Galactic Centre Deep Exposure has been 
analysed for the diffuse continuum from the Galactic 
ridge. A new catalogue of sources from the INTE- 
GRAL Imager IBIS has been used to account for 
their contribution to the celestial signal. Apparently 
diffuse emission is detected at a level ~ 10% of the 
total source flux. A comparison of the spectrum of 
diffuse emission with that from an analysis of IBIS 
data alone shows that they are consistent. The ques- 
tion of the contribution of unresolved sources to this 
ridge emission is still open. 

Keywords: gamma rays; X-rays; Galaxy; interstellar 
medium; cosmic rays; compact sources; INTEGRAL; 
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1. INTRODUCTION 



The inner Galactic ridge is known to be an intense 
source of continuum hard X- and soft 7-ray emis- 
sion. Th e hard X-ray emission was discovered in 1972 
(Bleac h et all Il972j) and has subsequently been ob- 
served from keV to MeV energies by ASCA, Ginga, 
RXTE, OSSE, COMPTEL and most recently by 
Chandra and XMM-Newton. 

While the physical process (e+e^ annihilation) pro- 
ducing the positron line and positronium contin- 
uum is clear, for of the remaining continuum it is 
not, a lthough nonthermal bremsstrahlung is most 
likely l|Dogiel et all 1200231) : the imphed photon lu- 
minosity in the continuum is a few 10'^* erg s~^. 
Its origin in a point-sou rce population seem s un- 
likely at keV energies (iTanaka et all Il999|) since 
there are no known candidate objects with the re- 
quired spatial and luminosity properties. In addi- 
tion, high-resolution imaging in X-rays with Chan- 



dra ijEbisawa et a ll2001tl shows a tru ly diffuse com- 
Donent. Recentlv iHands et al.l l|2004D have used the 
XMM-Newton Galactic Plane Survey to show that, 
in the 2-10 keV band, 80% of the ridge emission 
(l=19°-22°) is probably diffuse, and only 9% can be 
accounted for by Galactic sources (the rest being ex- 
tragalactic sources). 

Non-thermal 7-ray emission in the interstellar 
medium i mplies a very high luminosity in energetic 
particles l)Dogiel et all 12002^). Emissio n mecha- 
nis ms have b e en dis cussed bv lValinia et alJ i_2000a b) 
and iTanakal l|2002f l. and a solution to the prob- 
lem of energetics has been suggested by|Dogie^_et_al 



(2002b) but the issue is quite open. IValinia et al 



([2000^) propose a composite model with thermal and 
nonthermal components from electrons accelerated 
in supernovae or ambient interstellar turbulence. At 
MeV energies the origin of the emission is also uncer- 
tain l|Strong et alll200d) . The study of this emission 
is a key goal of the INTEGRAL mission. The high 
spectral resolution combined with its imaging capa- 
bilities lead to high expectations for investigating the 
nature of this enigmatic emission. 

Apart from its intrinsic interest, reliable modelling of 
the diffuse emission will be essential for the study of 
point sources in the inner Galaxy, since it contributes 
a large anisotropic background against which they 
must be observed. 

The most directly comparable results to ours 
(with respect to energy range) are from OSSE 
on the COMPTO N Gamma Ray Observatory 
iPurcell et alJ l)199(J) used coordinated observations 
of sources in the Galactic centre region from the 
SIGMA satellite to evaluate the source contribu- 
tion. They concluded that after correction for these 
sources, the OSSE flux was consistent w ith that mea- 
sured by O SSE at 1=25° and 339°. iKinzer et al] 
(|1999l I2OOH) give an OSSE spectrum ioi I = b = 0, 
but without independent imaging observations, so 
that the point-source contribution remains uncer- 
tain. 



2 



A preliminary analysis of diffus e continuum emissio n 
using SPI data has been given in lStrong et alJl)2003j) . 
using the first cycle of INTEGRAL GCDE observa- 
tions. It was concluded that the diffuse emission was 
detected at a level consistent with previous experi- 
ments, but the systematic errors were large, due in 
part to the uncertainty in the point-source contribu- 
tion. Only four of the strongest point sources were 
included in that analysis, since at least at higher 
energies (>100 keV) the source contribution is not 
thought to be critical. At lower energies an adequate 
7-ray source catalogue (e.g. from INT EGRAL) was 
not ye t available, and it was stated in IStrong et al.l 
(2003) that point-source contamination was likely. 
Note that the adopted method of analysis has the 
property that sources which are unaccounted for are 
mainly 'absorbed' into the fitted instrumental back- 
ground rather than affecting the diffuse signal. 

Me anwhile, a new analysis of IBIS data 
(Lcb run et all l2004 iTerrier et all MM) has 

shown that indeed a large fraction of the total 7-ray 
emission from the inner Galaxy is due to sources, 
at least up to 100 keV. This work has produced a 
source catalogue containing 91 sources, which can 
now be used as input to the SPI model fitting, 
giving a much more solid basis for the analysis. This 
exploits the complementarity of the instruments 
on INTEGRAL for the first time in the context of 
diffuse emission. 



2. DATA 

The INTEGRAL Core Program l|Winkler et all 
I2OO3I) includes the Galactic Centre Deep Exposure 
(GCDE) which maps the inner Galaxy (330° < I < 
30°, -20° < 6 < 20°) with a viewing time of about 
4Ms per year. The full region is covered in one 
GCDE cycle, and there are two cycles per year. Data 
from the first two cycles are used for the study re- 
ported here. 

We use data from the SPI (INTEGRAL Spectrome- 
ter) instrument; description s of the instrum e nt an d 
performance are given in IVedrenne etal\ (|2QQ3); 
lAttic et al. (2003). The energy range covered by SPI 
is 20 keV - 8 MeV, but here we restrict the analysis 
to energies up to 1 MeV; above these energies the 
statistics are small and the analysis is more difficult, 
so is reserved for future work. The data were pre- 
processed using the INTEGRAL Science Data Cen- 
tre (ISDC) Standard Analysis software (OSA) up to 
the level containing binned events, pointing and live- 
time information. GCDE data from 33 orbital revo- 
lutions from 47 - 123 were used. 2348 pointings were 
used, as shown in Fig 1; the total exposure livetime 
is 3.34 X 10^ s. The exposure per pointing is typ- 
ically 1800s, but can be as low as 300s in cases of 
high telemetry losses in the early part of the GCDE. 
The energy calibration is performed using instrumen- 
tal background lines with known energies; while this 
is a critical operation for line studies (where sub- 
keV accuracy is required), for continuum studies a 



standard calibration (~1 keV accuracy) is quite ad- 
equate. Various energy binnings were used, depend- 
ing on the available statistics as a function of energy. 
Only single-detector events are used here. 

The instrumental response is based on exten- 
sive Monte Ca rlo simulations and parameterization 
ijSturner et all |2003); this has been tested on the 
Crab in-flight calibration observations and shown to 
be reliable to better than 20% in absolute flux at the 
current state of the analysis. 

Since INTEGRAL data are dominated by instru- 
mental background, the analysis needs to have good 
background treatment methods. In the present work 
the background ratios between detectors are ob- 
tained by averaging the entire dataset over time; this 
has been found to give results which hardly differ 
from taking 'OFF' observations, and has the advan- 
tage of much higher statistics (few OFF data are 
available at present). A check on the s ystematics in 
this b ackground approach was given in lStrong et al.l 
l|200l . 



3. METHOD 



We use the program spidifjit which fits the data to 
a linear combination of input astronomical skymaps 
(e.g. HI and CO surveys, emission models) and point 
sources, together with background components. The 
fit is performed by maximum-likelihood with one pa- 
rameter per component and energy range. The back- 
ground is fitted per pointing using a template for 
the ratios between detectors from the average over 
all observations. Hence the time-dependence of the 
background is explicitly determined from the data 
themselves on the assumption of constant detector 
ratios. 

Since the distribution of the emission is unknown but 
certainly correlated with tracers of large-scale Galac- 
tic structure and the diffuse interstellar medium, we 
include line-of-sight integrated HI and CO surveys 
(which trace the atomic and molecular gas) as basis 
models. In addition, since the positronium emission 
is an important source of continuum below 511 keV 
and is believed to include in addition a component 
more concentrate d towards the Galact i c centre than 
the gas tracers (iKinzer et al.L 120011 [ vlilne et al.L 
120021 iJean et al.ll 200a'Knodlsede r etaL jgOO^. and 
noting that iKin zer et al. (2001) and Milne ct al] 
l|2002D found no significant difference between the 
line and continuum annihiliation angular distribu- 
tions, we include also a Gaussian with FWHM 10° 
centred on I — b — 0. Note that each fitted compo- 
nent has an independent spectrum so that it can be 
separated in the fitting; the sum of the diffuse com- 
ponents then gives the total diffuse emission spec- 
trum of the inner Galaxy. Additionally we remark 
that a coded-mask telescope is only sensitive to flux 
contrasts, so isotropic emission is suppressed. 

For the sources we use the catalogue from the 
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IBIS analysi s desc ribed in iLebrun et alJ l)2004(l : 
iTerrier et alJ l)2004|) and Terrier (private communi- 
cation). This contains 91 sources; only the source 
positions are required since the source spectra are 
determined in the model fitting. 



4. RESULTS 



Fig 2 illustrates the temporal variation of the back- 
ground scaling factor as determined in the fitting, for 
the energy range 18-178 keV. Other energy ranges 
appear very similar. The increase between the 1st 
and 2nd GCDEs is clearly visible around pointing 
sequence 1600. High points represent increases due 
to end-of-orbit passes or other physical effects. Gen- 
erally the fluctuations are at the few percent level 
and show that the technique effectively determines 
the time-dependence of the background. Alternative 
background treat ments using predi ctions based on 
activation tracers l|.Tean et al.L l2003^ provide consis- 
tent background histories. 

Fig 3 shows the spectrum of total emission from the 
3 diffuse components, and the total from the sources 
included in the fitting. The diffuse emission is ap- 
prox imately 10% of the summed sources, cons istent 
with ILebrun et all l|2004D : iTerrier et"al] |2n04lK It is 
slightly lower than given in iStrong et all l)2003^ due 
to the additional sources included, but still within 
the error bars quoted in that paper. The positron- 
ium edge is detected, at low significance. (Note that 
the broad bins used here to study the continuum 
suppress the 511 keV line). No diffuse emission is 
detected above 500 keV, but this is still consistent 
within the error bars with the COMPTEL fluxes. 



For comparison we show in Fig 4 the scaled 
n^b==0) OSSE spectrum and results from RXTE 
iRevnivts cv, 20Q3J and COMPTEL (Strong e t all 
I1999D. The OSSE flu xes are scaled as described 
in IStrong et al.l l|2flfl3|) (continuum scaled by 0.5, 
positronium scaled by 0.3), since a precise compari- 
son is difficult because of the different instrument re- 
sponses. A scaling factor less than unity is expected 
between the l=b=0 OSSE flux per radian and the 
full inner radian flux measured by SPI because of 
the concentration of the emission towards the Galac- 
tic centre. The RXTE results for |6| > 2° are scaled 
to the SPI region using information provided by M. 
Revnivtsev (private communication). 



5. DISCUSSION 



compatible with or somewhat lower than found in 
previous work. The fact that we can now explic- 
itly include a large number of sources in the analysis 
makes it more robust. 

The question remains as to whether the apparently 
diffuse emission is really interstellar or is the su- 
perposition of sources below the present detection 
threshold. This issue can only be answered in the 
future with population synthesis models and study- 
ing the logN-logS of the detected sources. We note 
however that at energies just below those visible to 
INTEGRAL, similar issues are now being addressed 
by Chandra and XMM- Newton (see Introduction). 
The finding of Hands et al. ( 2004) with XMM that 
~80% of the Galactic 2-10 keV emission is diffus e 
supports the RXTE analysis of iRevnivtsevI l)2003(l . 
Hence we do expect diffuse emission at some level 
in the INTEGRAL range unless a dramatic cutoff 
occurs over a factor of a few in energy. 

The apparent difference in the diffuse-to-source ratio 
between X-ray and gamma-ray energies can proba- 
bly be understood as a sampling effect: while the 
present gamma-ray ratio includes all the bright inner 
Galaxy sources, the XMM ratio mentioned above is 
for a small region which does not sample the bright- 
est sou rces. Using RXTE sc an observations, the re- 
sult of 'Ebi sawa et alJ ijlQQQ*^ shows that for a wide 
latitude region the X-ray keV diffuse-to-source ratio 
is comparable with what we find for the inner Galaxy. 

Evaluation of the combined results from X-ray and 
gamma-ray instruments will be essential and will in- 
deed be possible in the near future. 

In Fig 5 we compare the summed source spectra with 
the MeV ridge emi ssion measured by COMPTEL 
l|Strong et al.L IT999I) . This shows that a significant 
contribution to the 1-30 MeV emission could be 
due to the extension of the source spectra into this 
range, consistent with the requirement for an unre- 
solved source component suggested bv IStrong et al.l 



6. COMPARISON WITH IBIS RESULTS 



Fig 6 compares the diffuse spectrum from this work 
with t hat derived u s ing IN TEGRAL /IBIS ISGRI 
data bv lTerrier et all l|2004D (see Introduction). The 
results are consistent, both in the low-energy range 
(20 - 40 keV) , where the IBIS analysis gives an ex- 
cess not accounted for by detected sources, and at 
higher energies (40 - 220 keV) where IBIS was only 
able to derive upper limits. 



Detailed quantitative comparison with previous mea- 
surement of the diffuse continuum emission (from 
OSSE, RXTE, COMPTEL) is not possible because 
of the different sky areas and energy ranges covered, 
and the very different instrumental characteristics. 
Hence the comparison shown in Fig 4 is illustra- 
tive only, but shows that the INTEGRAL results are 



7. CONCLUSIONS 



We have exploited the complementarity of the in- 
struments on INTEGRAL to improve the SPI anal- 
ysis of diffuse Galactic 7-ray emission from 20 -1000 
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INTEGRAL GCDE Pointings Rev 47 - 123 
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Figure 1. INTEGRAL pointings used for this analy- 



sis. 



keV. Inclusion of the sources found by IBIS improves 
the sensitivity and rehability of the analysis. Dif- 
fuse emission is detected at a level rather lower than 
previously when the IBIS sources are accounted for. 
However this marks only the beginning of what will 
be possible by jointly analysing the data from both 
INTEGRAL instruments with their unique combina- 
tion of characteristics. 
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Figure 2. Instrumental background scaling factor de- 
termined by the fitting procedure. Note that the point- 
ing number is sequential, and there is a time gap be- 
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Figure 4- Diffuse continuum from this SPI analysis 
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spectr um around I — b — from \Kinzer et al\ 
\199!\ ). with components: exponentially cutoff power 
law (red), high-energy continuum power law (blue), 
positronium (gre en) and total (Mack ). Components 
are scaled as in \Strona et a,l,\ l200,i) and desc ribed 
in the text. 2. RXTE from \Revnivtsei\ \2003i) and 
Revnivtsev (private communicatio n) (lig ht blue bars), 
3. COMPTEL from \Strona et al\ (magenta). 
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Figure 5. Summed SPI source spectra from this 
analysis (red squares) compared with ridge e mission 
measured by COMPTEL from .Strong et aL 1^199!^) 
(magenta). 
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Figure 6. Diffuse Ridge emission measured by SPI 
(this work): dark blue, c ompared with results from 
IBIS \Terrier et al\.\20nA) : red. 
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